RNA Biology

Bacterial small RNAs (sRNAs) regulate gene expression by base pairing with target mRNAs, yet transcriptome-wide interactome mapping has shown that many sRNA-mRNA interactions detected in vivo have modest or no regulatory effect using orthogonal reporter assays. The features that determine functional outcome remain poorly defined. Here, we integrated Hfq-CLASH interactome mapping with matched transcriptomic and proteomic profiling in Escherichia coli and developed an interpretable machine-learning framework to identify the determinants that distinguish functional from non-functional interactions. Using sequence, structural, thermodynamic, duplex and protein-occupancy features, transcriptomic and proteomic responses were predicted with above-chance performance, achieving AUCs of 0.78 and 0.74, respectively. Feature attribution revealed that physical pairing alone is insufficient for regulation; instead, regulatory outcome is shaped by a coordinated interplay between RNA secondary structure, thermodynamic accessibility and local protein-binding context. Target-side Hfq occupancy emerged as a positive predictor of functional regulation, whereas AR2-domain occupancy on the sRNA was associated with non-responsive interactions, suggesting that distinct ribonucleoprotein states may separate productive regulation from non-productive binding. These findings indicate that the regulatory fate of an sRNA-mRNA interaction is an emergent property of its biophysical context and protein-binding environment, rather than a direct consequence of physical pairing alone.

Transfer RNAs (tRNAs) are key molecules that deliver amino acids to the translating ribosome according to their cognate codon encoded in messenger RNA (mRNA). Due to the modified nature of tRNA nucleotides, accurate tRNA quantification can be tedious, especially when dealing with small sample inputs. Here, we took advantage of an RNA-dependent RNA amplification method using T7 polymerase to quantify tRNA abundance in low biological input. Our method, called T7 High-resolution original RNA (Thor)-tRNA-Seq, showed reproducible and quantitative measurement of low tRNA inputs. Thus, our Thor-tRNA-Seq is a robust and reliable approach for the quantification of tRNA in samples with small and precious biological material.

This study proposes a novel therapeutic strategy to suppress cancer growth by modulating the MDM2-p53 axis via Alternative Polyadenylation (APA). MDM2 normally promotes tumorigenesis by ubiquitinating and degrading the tumor suppressor p53. In cancer cells, preferential use of proximal polyadenylation signals (PAS) results in shortened 3'UTRs, allowing oncogenic transcripts like MDM2 to evade nuclear sequestration mediated by Inverted Alu (IRAlu) double-stranded RNA structures. We hypothesized that forcing distal PAS usage would elongate the MDM2 mRNA, promoting its nuclear retention and reducing protein translation, thereby restoring p53 activity. Using CRISPR-Cas9, we targeted and deleted the most frequent proximal PAS in the MDM2 3'UTR of A549 cells. Successful genome editing was confirmed via PCR. As expected, Western blot analysis showed a significant reduction in MDM2 expression in PAS-edited cells. However, experimental outcomes contradicted our initial hypothesis: edited cells exhibited higher viability under doxorubicin treatment compared to wild-type cells. Furthermore, despite decreased MDM2 levels, a concurrent reduction in phosphorylated p53 (p-p53) was observed. These unexpected results suggest that MDM2 3'UTR elongation may trigger a non-canonical regulatory mechanism that bypasses the traditional MDM2-p53 interaction. This study highlights the complexity of post-transcriptional regulation and suggests that APA-mediated gene modulation can induce unforeseen compensatory survival pathways in cancer cells, necessitating further investigation into the broader functional landscape of elongated 3'UTRs.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=79 SRC="FIGDIR/small/745202v1_ufig1.gif" ALT="Figure 1"> View larger version (21K): org.highwire.dtl.DTLVardef@52e3edorg.highwire.dtl.DTLVardef@1f21637org.highwire.dtl.DTLVardef@909ee8org.highwire.dtl.DTLVardef@b0c907_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOGRAPHICAL ABSTRACTC_FLOATNO C_FIG Crosslinking and immunoprecipitation followed by sequencing (CLIP-seq) is widely used to identify the RNA targets of RNA-binding proteins (RBPs). However, its application to non-canonical RBPs lacking canonical RNA-binding domains and frequently displaying low or transient RNA occupancy, is limited by low signal-to-noise ratios, high input requirements and error-prone ligation steps during library preparation. To overcome these limitations, we developed soniCLIP, a streamlined CLIP-seq workflow that replaces RNase-mediated RNA fragmentation with sonication and uses a ligation-free strategy for library construction. soniCLIP is optimized for reproducible identification of enriched RBP-associated RNA regions from limited starting material. We benchmarked soniCLIP against the widespread eCLIP approach and observed reproducible recovery of known RBP-associated regions and target recovery comparable to ENCODE eCLIP, while requiring only 10% (500 {micro}g) of protein input. We further applied soniCLIP to the glycolytic enzyme and non-canonical RBP pyruvate kinase M2 (PKM2). We identified 197 significantly enriched RNA regions and validated selected targets by RIP-qRT-PCR and in vitro binding assays. By combining reduced input requirements, high reproducibility, a shortened 3.5-day workflow and the elimination of gel-based purification, soniCLIP provides an efficient and robust approach for the identification of RNA targets of canonical and non-canonical RBPs.

Transfer RNA (tRNA) modifications are essential for translational accuracy and cellular adaptation to environmental changes. While short-term modification dynamics are well documented, the impact of prolonged stress exposure on the global tRNA landscape remains largely unexplored. This study provides the first systematic profiling of tRNA modifications in Saccharomyces cerevisiae following long-term exposure to distinct stress types: heat, suboptimal pH, oxidative stress (paraquat and diamide), osmotic stress (NaCl and KCl), and genotoxic stress (MMS). To this end, we used a broad-range UPLC-MS protocol to quantify global changes in tRNA modification and identify stress-specific signatures. The results revealed that long-term stress triggers a global reprogramming of the tRNA epitranscriptome in a stress-specific and time-dependent manner. Importantly, while our findings confirm the previously reported temperature-sensitivity of wobble uridine thiolation, we also identified a complete or partial loss of 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2U34) modification upon exposure to paraquat and pH stress. Furthermore, this loss of thiolation is accompanied by an accumulation of the non-thiolated 5-methoxycarbonylmethyl (mcm5U) precursor, suggesting a stress-dependent impairment of the thiolation pathway. Next, we applied MarathonRT-based tRNA sequencing and showed that these modification changes occur independently of tRNA isoacceptor abundance. To further integrate these results, we devised a modification deviation (MDm) index, which indicates that the observed reprogramming is primarily linked to events that are independent from changes in tRNA abundance. Together, this study provides a comprehensive atlas of tRNA modification dynamics under prolonged stress, addressing a critical gap in our understanding of RNA-based translational control and establishes the MDm index as a robust quantitative framework to decouple the influence of tRNA abundance from global modification signals, providing a necessary metric for the field to interpret epitranscriptomic reprogramming. TABLE OF CONTENTS GRAPHIC O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=108 SRC="FIGDIR/small/745200v1_ufig1.gif" ALT="Figure 1"> View larger version (41K): org.highwire.dtl.DTLVardef@2cf571org.highwire.dtl.DTLVardef@1a5f3c2org.highwire.dtl.DTLVardef@266d61org.highwire.dtl.DTLVardef@d7214f_HPS_FORMAT_FIGEXP M_FIG C_FIG

The maternal-to-zygotic transition (MZT) is a major developmental event during which inherited transcripts are remodeled and zygotic transcription is established. Although parentally inherited long noncoding RNAs (lncRNAs) are present in early embryos, they have been thought to be dispensable. We have identified more than 2000 inherited lncRNAs in zebrafish embryos, but how these RNAs participate in regulatory programs during early development has remained unexplored. Here, the inheritance of selected zebrafish lncRNAs spanning a broad expression range were confirmed at the pre-MZT stage and full-length sequences were captured by Direct RNA nanopore sequencing. We show that 30% inherited intergenic lncRNAs are preferentially associated with active enhancers, annotated as such in DANIO CODE, whereas non-inherited intergenic lncRNAs rarely overlap with enhancers. Perturbation of five inherited intergenic lncRNAs, individually, using antisense oligonucleotides reduced the expression of their respective neighboring genes at 2.5, 4.3, and/or 6 hours post fertilization, indicating that these RNAs act as positive local regulators during MZT. Together, these findings identify inherited intergenic lncRNAs as enhancer-associated regulators with elncRNA-like properties during early embryogenesis.

The majority of eukaryotes encode some intron-containing pre-tRNAs. Splicing of these pre-tRNAs requires a dedicated tRNA splicing machinery. The fungal and trypanosome tRNA ligase, Trl1, and the human RNA ligase, RTCB, catalyze an essential step in tRNA splicing. However, Trl1 and RTCB are nonhomologous and biochemically and structurally distinct from each other. Therefore, Trl1 could serve as a broad-spectrum antifungal and anti-trypanosomal target. While the functions and requirements of the three catalytic Trl1 domains have been extensively characterized in the model yeast Saccharomyces cerevisiae, the roles of Trl1 orthologs in pathogenic fungi remain unexplored. Here, we validate Trl1 as one of the few promising novel drug targets for the development of antifungal therapeutics. Functional analyses of the three Trl1 domains show that only the "sealing" domain is essential for growth and viability in Candida albicans and Aspergillus fumigatus. In contrast, the two "healing" domains are dispensable in these pathogenic fungi, suggesting the presence of redundant healing enzymes, unlike in S. cerevisiae. These findings indicate that only the sealing domain is a good drug target. Our analysis also shows that the Mucor enzyme, which only contains the sealing domain, is essential. Using a Caenorhabditis elegans infection model of C. albicans, we further demonstrated that inhibiting Trl1 expression protects worms during an established infection. In contrast to these fungal pathogens, we show that all three domains of Trl1 are essential in Trypanosoma brucei. Our findings show that the essentiality of the Trl1 sealing is conserved in important human pathogens and provides an impetus for future drug development. SIGNIFICANCEFungal infections are an important cause of human disease and death and difficult to treat and there is an urgent need to develop additional drugs. Based on studies in yeast, one promising target for antifungal drug development is the tRNA splicing pathway. Human tRNA ligase is fundamentally distinct from the fungal one. To investigate the possibility of developing tRNA ligase-targeting drugs, we investigated the function of the catalytic domains of fungal tRNA ligase in different fungal pathogens. Surprisingly, only the first domain is essential in these pathogens and yeast is not a good model fungus. In contrast, all three domains of Trypanosome tRNA ligase are essential. These findings provide an impetus for future drug development.

ABSTRACT: The bacterial small RNA (sRNA) OxyS is expressed in Escherichia coli during oxidative stress. The sRNA OxyS enhances cell survival by controlling genes involved in the regulation of hydrogen peroxide (H2O2) and iron-sulfur (Fe-S) cluster formation. Here, we used the MS2 affinity purification coupled with RNA sequencing (MAPS) technique to identify new target mRNAs of the sRNA OxyS. Our analysis revealed a significant enrichment of mepS mRNA, which encodes a peptidoglycan endopeptidase that promotes cell growth. Our results confirm a previous report on the sRNA OxyS repressing the translation of mepS. We also found that an {Delta}oxyS background facilitates the emergence of mutations, conferring increased resistance to the last-resort antibiotics polymyxin B and E (colistin), but only in the presence of the target mepS gene. This suggests that the translation repression of mepS by OxyS could prevent mutations in bacterial DNA during H2O2-induced oxidative stress. Moreover, we show that adding the antioxidant thiourea or sequestering iron in the {Delta}oxyS background effectively reduces the emergence of resistance against both polymyxin B and colistin. These results suggest that reactive oxygen species (ROS), in conjunction with intracellular iron, play a key role in driving the emergence of antibiotic resistance. Overall, our work underlines a mechanism of antimicrobial emergence implicating oxidative stress, intracellular Fe, and cell wall remodeling in E. coli. IMPORTANCE: This study uncovers an underexplored link between peptidoglycan remodeling and oxidative stress responses during exposure to antibiotics. By elucidating how MepS and the sRNA OxyS interact in the presence of polymyxins and oxidative stress, our study suggests that MepS may exert an anti-mutator function. The repression of mepS translation by OxyS seems to limit the emergence of antibiotic resistance driven by DNA mutations. Together, these findings suggest cell wall remodeling and oxidative stress response pathways as promising targets to enhance antibiotic efficacy and limit the emergence of resistance.

RNA-targeting CRISPR systems are commonly evaluated by RT-qPCR, but guide RNA binding can confound these measurements. We show that crRNA alone produces apparent knockdown without reducing target RNA abundance, whereas RNA sequencing remains unbiased and reveals guide-associated transcriptomic perturbations. A simple RNA denaturation step before reverse transcription restores accurate RT-qPCR quantification, providing practical guidance for RNA-targeting CRISPR analysis and guide design.

Very short RNAs play critical roles in modern biology, and are thought to have been crucial for genome replication during the origin of life. Next-generation sequencing is an essential tool for characterizing pools of small RNAs, but current library preparation methods suffer from strong size and sequence biases. Here we present tinyRNA-seq, an optimized library preparation method designed to minimize length- and sequence-dependent capture bias enabling the sequencing of RNA fragments as short as 2 nucleotides. We use degenerate adaptor regions to reduce ligation sequence bias and facilitate unique molecular identifier (UMI) installation. We benchmarked tinyRNA-seq against commercial kits using a model primordial RNA genome consisting of hundreds of defined oligonucleotides ranging from 2 to 12 nucleotides. tinyRNA-seq reproduced the input RNA distribution without the size and sequence bias of the commercial kits. tinyRNA-seq also enables the detection of de novo oligonucleotide generation, an important process for the origins of life. Applied to biologically derived small RNAs including miRNAs, piRNAs, and cityRNAs, tinyRNA-seq showed significantly lower capture bias and recovered a wider range of sequences than commercial kits. tinyRNA-seq may thus provide a more complete and quantitatively accurate representation of small RNAs from both biological and chemical sources. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=94 SRC="FIGDIR/small/743385v1_ufig1.gif" ALT="Figure 1"> View larger version (29K): org.highwire.dtl.DTLVardef@95ee64org.highwire.dtl.DTLVardef@155fb06org.highwire.dtl.DTLVardef@1d3665forg.highwire.dtl.DTLVardef@1e61404_HPS_FORMAT_FIGEXP M_FIG C_FIG

Circular RNAs (circRNAs) are an emerging class of RNAs with biomarker potential, but their detection in liquid biopsies is challenging due to low abundance. We developed Ouro-seq, a novel long-read sequencing protocol optimized for full-length circRNA recovery. Applied to urine, cervico-vaginal self-samples from cervical cancer patients, and plasma from lung cancer patients and controls, Ouro-seq recovered 2-5 times more and substantially longer circRNA molecules than conventional methods. Plasma contained predominantly exonic circRNAs, while urine and cervico-vaginal samples were dominated by previously undercharacterized intergenic circRNAs. We also identified extensive alternative circularization and splicing events. Functional analysis revealed distinct specialization patterns: exonic circRNAs showed enhanced miRNA sponging potential, while circRNAs from unplaced genomic scaffolds demonstrated greater peptide-coding capacity. This study establishes Ouro-seq as a valuable tool for comprehensive circRNA characterization in low-yield clinical samples and advances circRNA biology understanding with potential biomarker discovery and disease monitoring applications. MotivationWhile circular RNAs (circRNAs) constitute a minor fraction of total RNA, they may play critical roles in cancer development. CircRNA concentrations are typically too low for detection by Oxford Nanopore Long-Read Sequencing (LRS), particularly in samples with limited RNA content, such as liquid biopsies. Consequently, LRS-based circRNA analysis from liquid biopsies remains unexplored. To overcome these technical limitations, we developed an optimized circRNA enrichment method utilizing short-amplicon suppression, enabling circRNA profiling from urine, plasma, and cervico-vaginal samples.

Small RNAs (sRNAs) regulate the expression of many genes including those involved in antibiotic resistance and bacterial virulence, making them potential therapeutic targets. A molecule that binds an sRNA could interfere with its ability to bind its target mRNA and disrupt the regulation mechanism. Randomization and screening of natural arginine rich motif (ARM) peptides led to peptides capable of interfering with the sRNA MicF's ability to regulate ompF in Escherichia coli. Molecular dynamics simulations suggested that this effect was not a result of a direct disruption of the MicF-ompF interaction. Instead, the peptides interfere with binding of the chaperone Hfq, which is required for MicF-mediated regulation. Subsequent testing demonstrated peptide specificity for MicF over two other Hfq scaffolds and the ability to disrupt regulation of two additional MicF targets. Together, these findings support the use of synthetic ARMs as a potential tool for modulating sRNA function in bacteria.

In contrast to their important structural and regulatory functions, such as in the metabolism of cyanobacteria, genes encoding small proteins are often not well characterized. Cyanobacteria use redox equivalents and energy from oxygenic photosynthesis to produce organic carbon compounds from inorganic carbon (Ci) and organic nitrogen compounds from inorganic nitrogen sources. Therefore, the assimilation and metabolism of carbon and nitrogen are coordinated at multiple levels in cyanobacteria. Here, we analyzed the Synechocystis sp. PCC 6803 gene ssr3189 encoding a 55 amino acids protein. Orthologs were detected in 665 cyanobacterial genomes defining COG5794 in the Database of Clusters of Orthologous Genes. Homologs in several eukaryotic algae suggest that Ssr3189 is an important protein that originated in cyanobacteria, was retained in algae after endosymbiosis, but was lost in plants. Polynucleotide kinase assays validated Ssr3189 as an RNA-binding protein. Deletion of ssr3189 resulted in lower pigmentation, delayed growth, and alterations in the expression of genes encoding transporters for nitrogen and Ci, and metabolic enzymes. Metabolomic analysis revealed a substantial overaccumulation of glutamine and tricarboxylic acid cycle intermediates in the deletion mutant, and further differences in the amino acid and organic acid pools compared to the wild type. Co-immunoprecipitation analysis yielded ribosomal protein S21, enolase and the Cas6-1 endoribonuclease as the most strongly co-enriched proteins, together with all other ribosomal proteins and a small set of metabolic enzymes. These findings are consistent with observations that ssr3189 encodes the ribosome-associated protein cS24 and suggest that it connects translation with metabolic control, and, potentially, RNA decay. IMPACT STATEMENTDespite considerable progress in analyzing microbial genomes, there are still substantial numbers of uncharacterized gene functions. Here, we analyzed a mutant lacking gene ssr3189 that is widely conserved, but phenotypically uncharacterized in cyanobacteria. This gene is important for growth at the optimum temperature and essential at lower temperatures. In its absence, important metabolites were overaccumulated, while genes involved in nitrogen and Ci uptake were dysregulated. The encoded protein binds RNA and interacts with proteins involved in translation and metabolism. The findings are consistent with a function as a ribosomal protein bridging protein synthesis and the regulation of metabolism.

N2-methylguanosine (m2G) is widely found at multiple positions in tRNAs across the three domains of life. Tryptophan tRNA from Thermococcus kodakarensis contains m2G at position 67. We previously proposed that the tRNA m2G methyltransferase Trm14 is responsible for m2G67 formation in tRNATrp from T. kodakarensis, although Trm14 was originally identified as the enzyme catalyzing m2G6 formation in tRNACys in Methanocaldococcus jannaschii. Thus, it remained unclear whether Trm14 could also methylate G67. Here, we characterized archaeal Trm14. Biochemical analyses using recombinant T. kodakarensis Trm14 revealed that the enzyme catalyzes m2G formation at positions 6 and 67 in T. kodakarensis tRNACys and tRNATrp transcripts, respectively. Mass spectrometric analyses demonstrated the loss of m2G6 and m2G67 in native tRNACys and tRNATrp, respectively, from a T. kodakarensis trm14 gene disruptant strain, providing direct evidence for the dual-site specificity of T. kodakarensis Trm14. The growth phenotype of the trm14 gene disruptant strain was comparable to that of the wild-type strain. In contrast, a trm14/trm11 double disruptant, in which trm11 encodes the tRNA m2G10/m22G10 methyltransferase, exhibited severe growth retardation at 95 {degrees}C. This suggests that m2G6/m2G67 and m2G10/m22G10 cooperatively contribute to cellular fitness at high temperatures. Biochemical analyses revealed that Trm14 methylates all 46 T. kodakarensis tRNA transcripts. Furthermore, we found that recombinant M. jannaschii Trm14 methylated both positions. In contrast, the bacterial ortholog TrmN modified only position 6 in tRNA. Overall, this study expands our understanding of archaeal Trm14 by demonstrating its broader substrate specificity and the physiological significance of these modifications under hyperthermophilic conditions.

Intrinsically disordered regions are widespread in RNA-processing machines. In Escherichia coli, RNase E uses its intrinsically disordered C-terminal domain (CTD) to recruit RNAs to the N-terminal catalytic domain, including mRNAs targeted by regulatory small RNAs (sRNAs), but the basis of substrate recognition and specificity is unclear. We engineered a protease-cleavable RNase E and used split-CRAC to isolate RNAs crosslinked to the AR2 sub-domain of the intrinsically disordered CTD fragment. AR2 preferentially engaged mRNAs and was depleted of sRNAs and sRNA-containing hybrids, supporting recognition of the mRNA. AR2 contacts concentrated on accessible A-rich motifs surrounding ribosome-binding sites and start codons, and purified AR2 recognised this motif in vitro. AR2 also contacted an AUAA motif in the rne translation-initiation region, and AR2 deletion increased RNase E abundance implicating this interaction in autoregulation. These findings define a relatively short AR2-binding motif and are consistent with CTD interactions with the 30S subunit that may provide additional specificity for a subset of mRNA translation initiation regions. SIGNIFICANCE STATEMENTMost RNA turnover in bacterial cells is carried out by the RNA degradosome, yet how this molecular machine checks and selects RNAs for degradation remains incompletely understood. We show that an intrinsically disordered region of the degradosome enzyme RNase E, termed AR2, preferentially binds A-rich sequences near sites of translation initiation. Through recognition of this shared sequence feature in a common functional context, AR2 may help the degradosome recognise messenger RNAs as a functional class. AR2 also contributes to feedback control of RNase E expression by recognising its own messenger RNA.

The QxxR motif is evolutionarily conserved within DEAD-box RNA helicases, including Drosophila Me31B and human DDX6, which post-transcriptionally regulate gene expression during animal development. A pathogenic H372R substitution (QxHR to QxRR) in the QxxR motif of human DDX6 has been associated with various developmental defects, but how this motif contributes to DDX6-family protein function remains unclear. Here, we used Drosophila Me31B as an in vivo model to investigate the QxxR motifs developmental role. We generated a Drosophila strain carrying the corresponding H333R missense mutation in Me31B and characterized its effects on female fertility, embryonic viability, germline development, and Me31B-associated molecular pathways. The me31BH333R mutation reduced female fertility in a gene dose-dependent manner, with homozygous mutant females being sterile. Embryos from the mutant females also exhibited primordial germ cell defects. Despite these developmental phenotypes, the me31BH333R mutation did not significantly alter Me31B protein abundance, global ovarian transcriptome or proteome profiles, or representative germ plasm mRNA and protein localization. In contrast, bait-normalized IP-MS analysis revealed altered enrichment of selected Me31B-associated proteins, including increased association of known Me31B interactors Trailer hitch (Tral) and Ypsilon Schachtel (Yps). These findings establish Me31BH333R as an in vivo model for investigating the conserved QxxR motif and suggest that disruption of this motif compromises development not through broad changes in gene expression, but potentially through altered composition or regulation of Me31B-containing ribonucleoprotein complexes.

Background: RNA molecules form thermodynamic ensembles, but interpretation often requires a single representative structure. Existing base-pair centroid estimators assess agreement at the level of individual base pairs and do not directly target nesting depth along the sequence. Methods: We introduce Mountain Centroid, which minimizes expected squared mountain-profile distance, and derive dynamic programming algorithms with and without RNA pairing constraints. We also combine the Mountain Centroid objective with the base-pair centroid gain. Results: Across 21,254 RNAStrAlign sequences, Mountain Centroid had lower median normalized mean squared mountain distance (NMSMD) than minimum-free-energy (MFE) and base-pair centroid ({gamma} = 1) structures, whereas its median base-pair F1 was lower. Imposing RNA pairing constraints improved base-pair F1 for 59.35% of sequences and reduced it for 3.58%. At an illustrative weight, the combined objective had median base-pair F1 similar to MFE while retaining lower median NMSMD than MFE and all tested {gamma}-centroid settings. Conclusions: Mountain Centroid represents an RNA structural ensemble with a single secondary structure that reflects how nesting depth varies across nucleotide positions. Combining mountain-profile and individual-base-pair criteria allows their relative contributions to be varied.

Pseudouridines ({Psi}s) are highly enriched in conserved, functionally critical regions of spliceosomal snRNAs, particularly within the U2 branch site recognition region (BSRR), which contains six {Psi}s in humans and three in S. cerevisiae. To investigate how U2 BSRR pseudouridylation influences branch site sequence (BSS) recognition, we developed a large-scale, high-throughput screening system in S. cerevisiae that allows a library of pre-mRNAs with randomized BSSs to be spliced in distinct U2 BSRR pseudouridylation backgrounds. Screening and next-generation sequencing (NGS) revealed that different U2 pseudouridylation backgrounds exhibit distinct recognition patterns and efficiencies for specific BSSs. Notably, {Psi}42 and {Psi}44 generally enhanced splicing, whereas {Psi}38 alone, and in some contexts together with {Psi}35, consistently impaired BSS recognition. The differential effects were validated using endogenous S. cerevisiae genes. In addition, splicing assays with engineered pre-mRNA constructs guided by the screening results demonstrated that BSRR {Psi}s directly influence BSS selection, supporting a model in which U2 pseudouridylation modulates BSS recognition and could contribute to alternative splicing in more complex eukaryotes. Finally, synthetic-lethality analyses with a Prp5 mutant, together with Prp5-U2-pre-mRNA binding assays, indicate that U2 BSRR {Psi}s are critical for Prp5 recruitment and, consequently, for proper U2-BSS interactions during spliceosome assembly. Collectively, these findings establish U2 pseudouridylation as a key determinant of branch site recognition and spliceosome function.

The ribosome, long regarded as a passive, uniform machine, has only recently been recognised as a direct regulator of translation. Mass spectrometry and sequencing approaches have shown that heterogeneity in ribosome composition exists, which can actively regulate the translational process. One source of this heterogeneity is the modification of ribosomal RNA (rRNA), primarily pseudouridylation (pseU) and 2'-O-methylation (2OMe), mediated by specific H/ACA and C/D box small nucleolar RNAs (snoRNAs). Here, we investigate how the stoichiometry of rRNA modifications varies during embryonic stem cell differentiation. Using the modification basecalling capability of Nanopore direct RNA sequencing, we have identified distinct stoichiometric changes in modification patterns between pluripotent and differentiated cells, revealing highly dynamic, site-specific regulation. Further, profiling of snoRNA expression during trilineage differentiation revealed differential expression of H/ACA and C/D box snoRNAs responsible for a subset of these dynamic modifications. By integrating rRNA and snoRNA sequencing approaches, we have built a comprehensive profile of rRNA modification dynamics during early embryonic cell fate decisions, highlighting potential regulatory mechanisms for ribosome heterogeneity during development. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=191 HEIGHT=200 SRC="FIGDIR/small/743918v1_ufig1.gif" ALT="Figure 1"> View larger version (26K): org.highwire.dtl.DTLVardef@48ca91org.highwire.dtl.DTLVardef@eb0418org.highwire.dtl.DTLVardef@159fc8corg.highwire.dtl.DTLVardef@d34f19_HPS_FORMAT_FIGEXP M_FIG C_FIG

Mutant tRNA mistranslation is a phenomena in which specific tRNA gene mutations cause translational errors wherein amino acids incorporated during translation differ from those coded by mRNA sequences. The biological consequences of mutant tRNA mistranslation are highly complex and contextual, but are often deleterious. Importantly, many variables exist that are likely to affect the outcomes of any one mutant mistranslating tRNA; biochemical characteristics of the exchanged amino acids, frequency of use, identity of mistranslated products, and so on. Here, we generate a carefully curated array of tRNA mutants to assess which characteristics of mistranslating tRNA variants are most predictive of deleterious phenotypes. We find that toxicity often arises from mutant tRNAs that induce dramatic changes in biochemical properties between exchanged amino acids, as well as mistranslation that occurs more frequently. Importantly, exceptions are also observed to each of these general rules, implying instead that some deliriousness may arise from more granular, product-specific mechanisms. Additionally, quantification of mistranslation via mass spectrometry demonstrates a weak relationship between quantity of mistranslated products and severity of toxic phenotypes. Lastly, despite previous establishment of mutant mistranslating tRNA models in Drosophila melanogaster, incorporation of higher frequency mistranslating tRNA variants was largely unsuccessful, and incorporation of lower frequency mistranslating tRNA variants produces no detectable developmental phenotypes. These experiments support the notion that although general rules may be capable of reasonably predicting their consequences, each mistranslating tRNA variant warrants individual consideration and investigation for thorough understanding of its biological outcomes and precise mechanisms of toxicity.